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Abstract. We formulate a theory combining the principles of a scalar- 
tensor gravity and the Rastall proposal of a violation of the usual conser¬ 
vation laws. In the resulting Brans-Dicke-Rastall (BDR) theory the only 
exact, static, spherically symmetric solution is a Robinson-Bertotti type 
solution besides the trivial Schwarzschild one. The PPN constraints can 
be completely satisfied for some values of the free parameters.The cosmo¬ 
logical solutions display, among others, a decelerate-accelerate transition 
in the matter dominated phase. 

1. Introduction 

The Brans-Dicke theory, the paradigma of scalar-tensor theories, is considered as 
an important alternative to the theory of General Relativity (GR) [Brans&Dicke, 
1962]. In this theory the gravitational coupling G is considered as a dynamical 
quantity represented by the field cj) which is introduced in the gravitational action 
through a kinetic term and a non-minimal coupling with the usual Ricci scalar. A 
new parameter u) quantifies the interaction of the scalar field and the gravitational 
term, such that as w ^ oo the General Relativity theory is recovered. Recent 
estimates using the PLANGK data indicates a value oj ~ 1000 [Avilez&Skordis, 
2013]. In spite of those observational constraints, small - or even negative - 
values of the parameter cj may be very interesting. They arise, for example, in the 
string theories in their low-energy limit [Lidsey et al, 2000]. When negative values 
of OJ are allowed, primordial singularity-free solutions emerge naturally from the 
Brans-Dicke theory [Gurevich et al, 1973]. Late time accelerated solution can 
also be achieved, but at the price of a negative gravitational coupling [Batista et 
al, 2001]. 

Some generalisations of GR evoke the gravitational anomaly effect, viz. 
Rastall’s theory. [Rastall, 1972; 1976] These generalisations touch one of the 
cornerstones of gravity theories: the conservation laws encoded in the null di¬ 
vergence of the energy-momentum tensor. Since the concept of energy in GR is 
an object of discussion, the possibility that the energy-momentum tensor has a 
non-zero divergence should be considered in some situations. The idea of viola¬ 
tion of the conventional conservation laws in the context of Brans-Dicke theory 
has been considered by Smalley [Smalley, 1974]. In this approach, the Klein- 
Gordon type equation for the scalar field was kept as in the Brans-Dicke theory 
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while the Einstein equations were changed accordingly. Here, we would like to 
revisit this proposal following a different path: we try to write down the field 
equations in such a way that the Brans-Dicke, GR as well as the ordinary the 
Rastall’s theory are recovered: we keep the violation of the energy-momentum 
tensor in a spirit very close to the original formulation of the Rastall’s theory, and 
the Klein-Gordon type equation as well as the Einstein’s equations are modified 
accordingly. 

We investigate the resulting theory in two situations: spherically symmetri¬ 
cal and cosmological configuration. In the former case, we obtain that the only 
non-trivial solution is represented by the Robinson-Bertotti metric [Bertotti, 
1959; Robinson, 1959] (its interpretation, however, differs from the conven¬ 
tional one). A solution that represents a star-like configuration is the “trivial” 
Schwarzschild one. At cosmological level, we show that accelerated solutions 
are possible in the dust phase of the cosmic evolution without introducing dark 
energy. We display a particular case where a decelerated/accelerated transition 
in the recent universe is achieved with a positive effective gravitational coupling. 


2. The theory 


The main idea of Rastall’s theory [Rastall, 1972] is the assumption that in curved 
space-time the usual conservation laws used in GR are violated. Hence, there 
must be a connection between the divergence of the energy-momentum tensor 
and the curvature of the space-time. According to this program, the divergence 
of the energy-momentum tensor may be written as 


— - _ — 

“ levrG • 


( 1 ) 


In equation ([T]) A is a free parameter codifying the deviation from the conser¬ 
vation. When A = 1 the traditional conservation laws are recovered. Equation 
([T|) is a phenomenological way to implement the gravitational anomaly due to 
quantum effects. 

In the context of the Brans-Dicke theory, we can make the identification 
G oc Hence, 


~ IhTT 
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Let us generalize Rastall’s version of the field equations to the Brans-Dicke case. 
Following the original Rastall’s formulation in the context of GR, we write, 

Rfiu ~ ~g^j.uR = ~ • (3) 

Gombining the hypothesis m and jS]), and using the Bianchi’s identities, we 
obtain that the scalar field (p must obey the equation: 


np = 
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3A - 2(1 - 2A)a; ^ 3A - 2(1 - 2A)w p 
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Equations mm form our complete system in this new formulation. When 
A = 1, the usual Brans-Dicke theory is recovered. 

The effective gravitational coupling today reads: 

2[2A+j3A-^l 

3A-2(l-2A)a; f ^ ’ 

When A = 1 we obtain the corresponding expression for the Brans-Dicke theory. 


3. Spherically symmetric static vacuum solutions 


The classical tests of theory of gravity are based on the motion of test particles in 
the geometry of a spherically symmetric object like a star or a planet. Hence, to 
verify the viability of the theory proposed, it is crucial to look for a spherically 
symmetric solution. As a first step, the (exterior) solution representing the 
space-time of a star-like object is considered. 

In the vacuum case, the equations reduce to 


R'^ = 0 , 




+ 

□0 = 


U) 

T 2 


A 


2(1 - 2A) 
(1 + A) 


9tiu4’\p4* 


IP 


+ 


' 2(1-2A) 
a;(l - A) (t>'P(f).p 
'3A - 2(1 - 2A)a; cj) 




( 6 ) 

(7) 

( 8 ) 


The first of these equations leads to i? = iio = constant. Hence, in vacuum 
the Ricci scalar is necessarily constant. The case Rq = 0 corresponds to the 
Schwarzschild solution of GR. 

Let us consider a metric in the form: 


^ g27^^2 _ _ g 2 / 3^^^2 ^ ^-^2 


(9) 


The functions a, /3 and 7 depend on the radial coordinate r only. First we find 
that the constant Rq is given by: 


Rn= uj 


3-h2a; 


3A-2(1 -2A)w 
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Now, the D’Alambertian reads: 


□ (/) = 




V^9 


= —e 


^—2arin 


[</." + (7' + 2/l'-a')0'], 


( 10 ) 


( 11 ) 


where g = 

Let us choose the radial coordinate such that a = 7 -|- 2/3. Writting the 
Einsteinian equations as 


- ^2 
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we obtain, in the extended form: 

7 + -77 = 


-uj- 
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7 " + 2/3" - 2/3'(/3' + 27 ') - ^( 7 ' + 2/3') = 

/3" + /3'^_e2(7+/3) = 
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Equations (I8I10I) lead to two supplementary equations: 

3 + 2w 


Ro = -oj 


3A - 2(1 - 2A)a; 
1 — A (t 


e I — I 


/ \ 2 


(kJ ’ 


/2 


= -UJ- 
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The only self-consistent solution for the above equations are: 

a = ao-ln(r/ro), 

7 = 7 o-ln(r/ro), 


= /3o = ^(ao- 70 ), 

1 

0 = (f>o(r/ro)i-^, A = -ca 


1-A 
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Hence, the metric is 
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ds^ = e 




(r/ro)^ 


— e 


270 . 




{r/rof 


— e 


0:0—70 




(16) 

(17) 


(18) 

(19) 

( 20 ) 
( 21 ) 


( 22 ) 


If the scale rg is chosen such that rg = e“° making redefinitions t 
e-(70+ao)/2^ s ^ e-iio-ao)/2^ ^ —7> rro, we arrive at: 


ds^ = 



(23) 


which is the so-called Robinson-Bertotti solution that is obtained, in the context 
of GR, by considering an electromagnetic field. Hence, no black hole solution 
is possible. This solution appears as the only non-trivial (non-Schwarzschild 
solution) vacuum solution. 

A PPN analysis reveal that the classical tests of a gravitation theory are 
equally satisfied as in Genera Relativity if A = 0 [Garames et al, 2014]. This 
makes this theory quite competitive. 


4. Cosmology 

Let us consider an isotropic and homogeneous space-time described by the flat 
Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, 

ds^ = dt^ — a{t)^{dx^ -|- dy^ -\- dz^)^ 


(24) 















5 


and an equation of state of the type p = ap, with a = constant. In this case the 
equations of motion read: 


3 
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Equations (I25l) - (l28p form a rich and complex system. In order to get a 
hint on which kind of solutions they predict, we consider power-law solutions, in 
the first place. The power-law solutions constitute a very restrictive case, but 
they can indicate the kind of cosmological solution we can expect from the BDR 
theory. Hence, suppose the solutions have the form a = aot^,4> = P = 
where oq, y>o, Po, s, p and q are constants. This system admits eight pairs of 
roots for {s,p). For the dust case, a = 0, one of the pairs corresponds to the 
Minkowski case, p = s = 0. Another one is s = p = 1/2. A third root implies 
a curious configuration with s = 0 and p = 2, that is, a static universe, with a 
varying gravitational coupling. Among the other five pairs, two incorporate an 
accelerated regime of the expansion while remaining three describe a decelerating 
universe. The overall situation is described in reference [Caramfe et al, 2014], 
where the meaning of these different results are analysed. 

The field equations depends not only on the values of A and ui, but also 
on the value of the initial conditions. We look for an example of a decelera¬ 
tion/acceleration transition during the matter dominated phase (a = 0). Figure 
1 shows the behaviour of the Hubble function H = ^ and deceleration parameter 

q = —1 — jp, for oj = 1 and A = — 1, undergoing this transition. Note that the 
effective G > 0 stays positive (see ©)• 


5. Conclusions 

In this work we have combined the idea of a scalar-tensor theory of the Brans- 
Dicke type and Rastall’s proposal of a gravitational anomaly encoded in the 
violation of the conventional conservation law for the energy-momentum tensor. 
In doing so, we end up with two free parameters: the usual Brans-Dicke param¬ 
eter OJ and Rastall’s parameter A, representing a degree of the non-conservation. 
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Figure 1. The left panel represents the first root for s displaying 
acceleration of the scale factor, the red line separating positive (left) 
and negative (right) values. The left center panel represents the first 
root for p, the red line separating positive (second and fourth quadrants) 
and negative values (first and third quadrants). The behaviour of the 
Hubble function H and deceleration parameter q for w = 1, A = — 1 
and a = 0 are displayed in the right center and right panels. 


We have investigated the BDR theory in two contexts: spherically symmetric 
static solutions and cosmological regime. In the first case, we found that the only 
possible non-trivial analytical solution is a Robinson-Bertotti type solution. The 
only possible solution in the BDR theory that can represent a star is the usual 
Schwarzschild solution corresponding to the trivial configuration where the scalar 
field is constant. 

For the cosmological case, we found power law solutions for the matter 
dominated phase, some of them representing an accelerating expansion, others, 
decelerating. We have shown that a decelerating/accelerating transition can be 
achieved in the matter dominated phase in the BDR theory. 
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